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ABSTRACT 

We present results from an optical imaging and spectroscopic monitoring 
campaign on the young, low-mass eruptive variable star V1647 Orionis. The 
star and associated nebulosity (McNeil's Nebula) were observed over the period 
February 2004 to February 2006 with observations commencing a few months 
after the original outburst event occurred. Using the Gemini North telescope, we 
obtained multi-band optical imaging photometry and medium-resolution long- 
slit spectroscopy of V1647 Ori on an approximately monthly interval. During 
this period, V1647 Ori remained at, or close to, peak brightness and then faded 
by 5 magnitudes to close to its pre-outburst brightness. This implies an outburst 
timescale of around 27 months. Spectral features seen in both emission and 
absorption varied considerably during the monitoring period. For example, the 
Ha line changed significantly in both intensity and profile. We present and 
discuss the observed photometric and spectroscopic changes and consider how 
this eruptive event relates to the early formative stages of low-mass stars. 

Subject headings: stars: individual(V1647 Orionis) - Reflection nebulae - Accre- 
tion, Accretion disks 



1. INTRODUCTION 



When the amateur astronomer Jay McNeil discovered a new nebula in the L1630 molec- 
ular cloud in Orion in January 2004 (McNeil 2004), little did he know that it would spark 
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significant and extensive worldwide follow-up investigations. The nebula he discovered has 
subsequently been designated "McNeil's Nebula" and its illuminating/exciting star named 
V1647 Orionis (Samus 2004). To date, around 30 publications have resulted from these 
studies spanning the electromagnetic spectrum from X-ray to radio wavelengths. In these 
research papers, considerable speculation has been put forth as to the nature of this erup- 
tive event, yet it is generally agreed that it is the result of a rapid and massive increase 
in accretion onto the surface of the young star. The reader is referred to the papers by 
Abraham et al. (2004a,b, 2006), Acosta-Pulido et al. (2007), Andrews, Rothberg, & Simon 

(2004) , Aspin et al. (2006), Aspin, Beck, & Reipurth (2008, herein ABR08), Aspin, Greene, 
& Reipurth (2009), Aspin et al. (2009), Briceno et al. (2004), Brittain et al. (2007), Fedele 
et al. (2007a,b), Gibb et al. (2006), Grosso et al. (2005), Kastner et al. (2004, 2006), Kospal 
et al. (2005), Kun (2008), McGehee et al. (2004), Mosoni et al. (2005), Muzerolle et al. 

(2005) , Ojha et al. (2004, 2006), Reipurth & Aspin (2004), Rettig et al. (2005), Semkov 
(2004, 2006), Tsukagoshi et al. (2005), Vacca, Cushing, & Simon (2004), Vig et al. (2006), 
and Walter et al. (2004) for further discussions. 

Such eruptive events have been observed in a number of young stars dating back to the 
observations and discussions by Herbig (1966). In this and subsequent papers, Herbig (1966, 
1977, 1989) suggested that these eruptions can be classed as either short-term EXors events 
(after the progenitor EX Lupi) lasting between a few weeks to a few years, or long-term 
FUors events (after the progenitor FU Ori) lasting decades to possibly even a century. Much 
emphasis has been placed on determining what type of outburst has occurred in the case 
of V1647 Ori, yet more compelling perhaps is understanding whether the FUor and EXor 
classes are fundamentally similar in their origins or whether they are the result of distinct 
phenomena occurring on different timescales. 

If the outburst of V1647 Ori is a FUor type eruption, then it is the first in almost 40 
years since that of V1057 Cyg (Welin 1971). A large-scale effort to monitor V1647 Ori was 
therefore undertaken involving use of both Gemini Observatory 8-meter telescopes (North 
and South) and the NASA IRTF 3-meter telescope over three observing semesters and 
utilizing six different facility instruments. This campaign resulted in approximately monthly 
optical imaging and spectroscopic and near-IR spectroscopic observations of V1647 Ori. In 
this paper, the first in a series describing these data, we present optical observations spanning 
the outburst phase, from the first follow-up observations taken soon after its discovery (taken 
in 2004 February and published in Reipurth & Aspin 2004) to the time when the source had 
faded to its pre-outburst brightness (in February 2006). Clearly, the approximately two year 
lifetime of this event suggests that, at least superficially, the eruption of V1647 Ori is more 
similar to EXor event than those occurring in FUors. 
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In this paper we present optical imaging and spectroscopic observations of V1647 Ori 
taken between February 2004 and February 2006 and discuss the variations present in these 
data. The observations and data reduction are described in §2. In §3 we present details 
of the temporal changes observed, and in §4 we consider what these observations imply 
regarding two physical components of the outburst, i.e. the fast wind and the accretion 
process. Finally, in §5 we piece together a time-line for the event and detail the evolution of 
the eruption as implied by the data from our monitoring campaign. 

2. OBSERVATIONS & DATA REDUCTION 

All but one dataset presented below were acquired using the "Frederick C. Gillett" 
Gemini North telescope located on Mauna Kea, Hawaii using the facility optical imager and 
spectrograph, GMOS-N (Davies et al. 1997; Hook et al. 2004). The imaging observations 
used standard GMOS-N g', r', i', and z' filters which were designed to be as close as possible 
in characteristics to the Sloan Digital Sky Survey (SDSS) filters (Fukugita et al. 1996). 
For the spectroscopic observations, all but one used the blue 600 lines/mm grating and 0'.'5 
wide long-slit. The one exception, the first observation taken, used the red 831 lines/mm 
grating, again with a 0'.'5 wide long-slit. The spectral resolution of data taken were therefore 
~4400 and ~1700, corresponding to 0.34 and 0.45 A pixel -1 respectively. The complete 
observation log is presented in Table [TJ In total, we acquired photometry at 16 different 
epochs and spectroscopy at 15 different epochs over the eruption period. 

The additional spectrum of V1647 Ori was acquired using the W.M. Keck II telescope 
located on Mauna Kea, Hawaii, using the facility high-resolution echelle spectrograph HIRES 
(Vogt et al. 1994). The observations were made on UT 2004 September 24 when the star 
had an optical V magnitude of about 17. The spectra covered the wavelength range 5800 to 
7150 A and 7320 to 8700 A at a nominal spectral resolution of R~46,000 using a 0'.'86 wide 
slit. The on-source exposure time used was 60 minutes. 

2.1. Imaging data reduction and calibration 

All datasets were reduced using the Gemini IRAF data reduction package vl.8. Specifi- 
cally, we used the routines gireduce for basic instrument signature removal (trim the images, 
subtract the master bias image, and divide by the normalized master twilight flat field), and 
gmosaic for combining the data from the three GMOS-N CCDs into one image. In addition, 
multiple exposures of the region were coadded using the IRAF routine imcombine. Aperture 
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photometry of V1647 Ori was performed using the Starlink program gaia (Draper et al. 
2008). An aperture diameter of 2" was used for the photometric calculation and, due to 
the bright nebulosity immediately surrounding the star, the sky signal was estimated from 
numerous apertures located in blank sky regions throughout the image. The full width half 
maximum (FHWM) seeing values in the acquired images are also given in Table [TJ 

Since some of the data were acquired during non-photometric conditions, a photomet- 
ric calibration was achieved relative to a dataset taken on a photometric night, UT 2004 
February 14. A series of six field stars were used to provide a calibration sequence; their 
position and magnitudes are listed in Table [2j Since the region containing V1647 Ori is one 
with active star formation and young stars are well-known to be variable, using six field 
stars as calibrators was considered an acceptable way of minimizing errors introduced by 
intrinsic variability of the calibrator stars themselves. These stars, together with V1647 Ori, 
are identified in Fig. [TJ Using instrumental magnitudes derived from calibrators, we ad- 
ditionally studied the variability of individual calibrators themselves and concluded that 
uncertainties resulting from source variability were small and within the quoted uncertain- 
ties on the V1647 Ori photometry. In addition to this, and to provide a consistency check, we 
have used the SDSS photometry of V1647 Ori from UT Nov 11, 1998. These data were part 
of the early " Orion" release (Finkbeiner et al. 2004) downloaded from Princeton Universit}0, 
and quoted in McGehee et al. (2004). Photometry of the same calibration sequence stars 
was extracted from the SDSS data and found to be consistent with the photometry from our 
GMOS images within the associated errors (shown in Table [2]). 



2.2. Spectroscopic data reduction 

Again, all the Gemini/GMOS data were reduced using the Gemini IRAF data reduc- 
tion package vl.8. Specifically, we used the routines gsreduce, gstransf orm, gsskysub, 
gscrrej, and gsextract for basic instrument signature removal (trim the images, subtract 
the master bias image, and divide by the normalized flat field), wavelength calibration, sky 
line subtraction, cosmic ray removal, and point-source spectrum extraction, respectively. 

The resultant spectra were not flux calibrated nor were atmospheric features removed 
since we were primarily interested in the spectral features present and their intrinsic variabil- 
ity. In addition, Ha has been shown to be a relatively insensitive diagnostic for accretion rate 
due to it, and other Balmer lines, becoming optically thick at high accretion rates (Muzerolle 
et al. 1998), thus negating any requirement for flux determinations. Fig. 8 of Muzerolle et 



see URL http://photo.astro.princeton.edu/oriondatarelease 
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al. (1998) show model predictions for Ha flux vs. accretion rate and it is clear that the 
relationship saturates at around 10~ 8 M yr _1 . Such rather low accretion rates are typically 
found in 'weak-line' T Tauri stars (WTTS) and low activity 'classical' T Tauri stars (CTTS) 
like DN and DQ Tau and not eruptive variables such as V1647 Ori. 

For our HIRES data, a standard reduction was performed, including bias correction, 
flat-fielding, scattered-light correction, order extraction, and wavelength calibrations using 
standard routines in the IRAF echelle package. 

3. RESULTS 

3.1. Optical photometry 

In Fig. [2] we show our optical photometry of V1647 Ori from 2004 February to 2006 
February for all four passbands, namely g', r', i', and z'. The horizontal dashed and dotted 
lines represent the pre-outburst SDSS brightness of the star (labeled filterpo). The solid 
lines simply join the discrete observation points as a guide to the eye. The two long 'gaps' in 
the observation sampling are the 2004 and 2005 summer months when Orion was not visible 
during the nighttime period. 

Over the summer of 2004, the source appeared to remain at approximately the same 
brightness, however, it is likely that some small-scale variability occurred. Our data from the 
2004-2005 winter months suggests that variability at the ~0.5 magnitude level was present. 
Such fluctuations were also reported by McGehee et al. (2004), Semkov (2004), Walter et 
al. (2004), Ojha et al. (2006), Acosta-Pulido et al. (2007), and Fedele et al. (2007) during 
several different observing periods. 

Over the summer of 2005, V1647 Ori exhibited a small yet significant decline in bright- 
ness, approximately the same in all filters, and amounting to <1 magnitude. Immediately 
after this period, in 2005 September, a major decline phase began. From 2005 September 
onwards, the trend was for monotonic dimming, with possibly some slowing occurring in 
late 2005. By the time of our last g' band photometric observation (in early 2006 January), 
V1647 Ori had already faded by about 4.4 magnitudes from its maximum brightness. At r', 
we obtained one further observation (in 2006 mid-February) which shows a total fading of 
over 5 magnitudes. At this time V1647 Ori was already close to its pre-outburst SDSS r' 
magnitude of 23.04. 

In late 2006 December we took a short R-band exposure of V1647 Ori at the University 
of Hawaii 2.2m telescope on Mauna Kea, Hawaii which showed that the source had remained 
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very faint and close to the above pre-outburst SDSS brightness. In addition, ABR08 pre- 
sented photometry from 2007 February which showed V1647 Ori still close to its pre-outburst 
brightness (r'=23.26±0.15). 

3.1.1. The GMOS and SDSS photometric systems 

Below, we make use of SDSS photometry and therefore we briefly consider the signifi- 
cance of any differences in photometric systems involved to give confidence to the analysis 
that follows. As we noted above, the GMOS g',r',i',z' filters were designed to be close 
to identical to those used by the SDSS survey (Smith et al. 2002). We therefore expect 
the photometric results from GMOS to be quite similar to those from SDSS. However, to 
quantify this we have studied the GMOS and SDSS photometry of the calibration sequence 
used to boot-strap the photometry of V1647 Ori over the monitoring period. In Fig. [3] we 
show the photometry of one representative field star labeled S4 in Fig. [TJ These photo- 
metric values were derived using standard GMOS zeropoints for the observing period (i.e. 
ZP(g')=27.93, ZP(r')=28.18, ZP(i')=27.90, ZP(z')=26.77, Jorgensen 2009) determined us- 
ing observations of Landolt (1992) standard stars observed close in time to the UT 2004 
February 14, V1647 Ori observations. The magnitudes of the Landolt standards (in the 
Johnson-Kron-Cousins photometric system) were transformed to the SDSS photometric sys- 
tem using the relationships derived by Smith et al. (2002). The horizontal lines through the 
observations are the mean magnitudes of S4 in each GMOS filter (dot-dashed lines) and the 
SDSS magnitude of S4 from November 1998 SDSS observations (dashed lines). We note that 
the slight shift between these lines suggest small differences in stellar magnitudes, however, 
the shifts are in all cases smaller than the associated la error on the data. We conclude, 
therefore, that the differences between the GMOS and SDSS photometric systems are not 
significant for the analysis presented below. 

3.1.2. Optical colors 

During the whole outburst period, when the source was at or close to its maximum 
brightness, the optical colors remained approximately constant. In Fig. [4] we show the trend 
in these colors over the monitoring period. The available SDSS pre-outburst colors (r'-i' and 
i'-z') are also shown as horizontal lines at the right edge of the plot. We note that during the 
outburst, i) the optical color indices are larger at shorter wavelengths, ii) both the r'-i' and 
i'-z' colors are bluer during outburst (by 0.4 and 0.5 magnitudes, respectively), and Hi) there 
is a clear trend for these colors returning to their pre-outburst values as the decline phase 
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progresses. In Fig. [5] we show the optical r'-i' and i'-z' colors plotted in a two color diagram. 
The colors of V1647 Ori after outburst, but before the decline had started, are the cluster of 
points near r'-i'~1.8 and i'-z'~1.5. The point labeled SDSS 1135 is from Modified Julian Date 
(MJD 2450000+) 1135 which is UT 1998 November 17 (the pre-outburst SDSS source color). 
The three points labeled 3702, 3729, and 3740 are the source colors at those MJDs and are 
from the decline phase. These points are significantly closer to the pre-outburst color than 
those during the outburst. In this plot, we also show the locus of main sequence dwarf colors 
(solid line) adapted from Fig. 1 of Finlator et al. (2000) and reddening vectors (dashed lines) 
for a ratio of total to selective absorption R=Ay/E(B-V)=3.1 extending from the extremities 
of the dwarf locus. These vectors represent a visual extinction Ay=5 magnitudes. We have 
used the tabular data of D. Finkbeiner (private communication]^ to calculate the effective 
change in r'-i' and i'-z' colors over the Ay range plotted. All the colors of V1647 Ori lie 
between the two reddening vectors. Dereddening the SDSS 1135 color into the dwarf locus 
suggests a spectral type of late-K to early M dwarf. During the outburst phase and through 
the decline phase, the colors of V1647 Ori first become more blue and then more red. This 
can be explained by either one, or a combination, of effects. Either a change in intrinsic 
source color or a color and extinction change has occurred. We explore this effect in a plot 
of reddening invariant colors vs. time shown in Fig. [6] Reddening invariant colors are, as 
the name suggests, colors that do not change as extinction along the line-of-sight changes. 
These have been used by McGehee et al. (2004) in their study of V1647 Ori using SDSS 
photometry and McGehee et al. (2004) in a study of accretion in low-mass young stars. 
Following Table 4 of McGehee et al. (2004), the reddening invariant color we use, Q r j 2 is 
defined as 

Q riz = ( r ' - i') - 0.987(?' - z) for R v = 3.1 (1) 

Here we chose the standard interstellar value of Ry=3.1, however, we note that the 
change in Q r ; 2 for somewhat larger grains, e.g. Ry=5.5, is relatively small (~10%) when 
compared to the associated uncertainties on the measurements. In Fig. [6] we perhaps see 
a slight trend in Q riz suggesting a larger value as the eruption proceeds into the decline 
phase. However, the change in value is a small effect with respect to the associated errors. 
If this difference is real, it would suggest that the intrinsic colors of the source are different 
out of outburst than in outburst. Such a color change is anticipated due to the added 
optical continuum luminosity from the enhanced accretion during the outburst (McGehee 
et al. 2004). Fig. 2 of McGehee et al. (2005) plotted numerical model predictions for the 



2 http:/ /www. astro. princeton.edu/~dfink/sdssfilters/ 
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Q r i Z of low-mass young stars vs. effective temperature, T e jf, and as a function of surface 
gravity, log(g). McGehee (private communication) has subsequently produced a version of 
this plot for us showing the variation in Q riz up to T e ff= 10000 K. We do not wish to 
over-interpret the current dataset due to the associated errors, however, comparing the Q riz 
value with the aforementioned plot, suggests that in its eruptive state the T e ff of the optical 
emission from V1647 Ori determined from its observed Q r j 2 value (~0.35±0.15) is in the 
range 3000-4000 K. Distinguishing between different surface gravities in this plot would be 
impossible since the range of Q riz values for log(g) of 3.5 and 5.5 at this T e ff is 0.25 to 0.35, 
respectively. Unfortunately, we note that the SDSS pre-eruption value of Q r ; 2 from McGehee 
et al. (2004)1 is of relatively poor quality (Q r j 2 =0.25±0.23) due to the faintness of V1647 Ori 
in r' in November 1998 (MJD 1135, r'~23.04±0.22) and is therefore not considered further. 



3.1.3. The pre-outburst to post-outburst light curve 

In order to produce as complete a light curve of the 2003-2006 outburst of VI 647 Ori 
as possible, we have combined our photometry with selected published results. In Fig. [7] we 
show our SDSS i' band GMOS photometry from Table [3] together with the lc photometry 
of Briceno et al. (2004). We have used the polynomial transformation equation of Ivezic et 
al. (2007) to estimate SDSS photometric i' values from the Kron-Cousins (Landolt 1983) lc 
photometry of Briceno et al. (2004). The transformation equations used was 

i' = I c + 0.0307(r - i'f - 0.1163(r' - i') 2 + 0.3341(r' - i') + 0.3584 (2) 

Since this (and any) transformation is source color dependent, we assume an intrinsic 
color for V1647 Ori of r'-i'=1.9, the mean of the pre-outburst and outburst colors. This 
results in a transformation correction of 0.79 magnitudes. Additionally, since the software 
aperture used by Briceno et al. is considerably larger than ours (aperture radii of 4'. / l and 
l'.'O, respectively), we have estimated from our UT 2004 February 14 image an aperture 
correction for the transformed Briceno et al. data. Since the nebulosity local to V1647 Ori 
was reasonably bright in Briceno et al. images from 2004 December 15 onwards, we apply 
the aperture correction to only their data points subsequent to this date. The correction 
used was 0.62 magnitudes. The result of the transformation and aperture correction is to 
give a reasonably good match (Am=0.1 m ) between our photometry (from MJD 3049) and 



3 Tablc 5 from McGehee et al. (2004) quotes an eruptive value for Q r i Z for V1647 Ori of 0.81. This used 
preliminary photometry from Reipurth & Aspin (2004). After a more thorough study of the photometry 
presented in that paper, the true eruptive phase value for Q r i Z is 0.35. 
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the closest value from Briceno et al. (from MJD 3036). Also shown in this figure is the lc 
photometry from Acosta-Pulido et al. (2007). We have used the same transformation to 
SDSS photometry given by Ivezic et al. (2007) and shown in Eqn. 2 above. In Fig. [7] our 
data are shown as open circles, the Briceno et al. data as filled dots, and the Acosta-Pulido 
et al. data as filled stars. We additionally note that the pre-outburst (MJD 1135) SDSS 
i' photometry from McGehee et al. (2004) (indicated by the black dotted line in Fig. [7]) is 
consistent with our last data point (MJD 3782) suggesting that V1647 Ori had more or less 
returned to its pre-outburst brightness in 2006 February. 

Another dataset that has good temporal coverage using only one telescope/instrument 
combination is the R band photometry of Ojha et al. (2006). They present 13 data points 
covering the period from peak brightness to approximately half way down the steep decline 
phase. These data are shown in relation to our r' band photometry in Fig. [8J Due to the 
extensive overlap of data, we have made no attempt to match the two photometric systems 
and have merely shifted the Ojha et al. data (filled dots) by +1 magnitude (fainter) to match 
our data (open circles) as well as possible. We again note that our final r' photometric point 
(MJD 3782) is very close in value to the pre-outburst (MJD 1135) SDSS value of 23.04±0.22 
from McGehee et al. (2004). 

We have additionally compiled all published near-IR K-band (2 ^m) photometry and 
plotted it together with the above r' light curve in Fig. [91 This comparison shows that in the 
near-IR the change in brightness of V1647 Ori was less than in the optical and amounted to 
~2.8 magnitudes from the October 1998 2MASS photometry (horizontal dashed line). The 
data were taken from Acosta-Pulido et al. (2007) and Ojha et al. (2006) although some of 
the data points in Ojha et al. (2006) were taken from other sources (specifically Reipurth 
& Aspin 2004, McGehee et al. 2004, and Ojha et al. 2005). Unfortunately, there is no 
K-band photometry from the spring and early fall of 2006 and therefore we cannot, with 
certainty, say whether the two light-curves are in phase or if there is a relative delay between 
them. The few points in late 2006 and early 2007 suggest that the major decline phase of 
V1647 Ori was somewhat shallower at K than in the optical although this relies on so few 
points that it is by no means certain. Nonetheless, the main fact arising from this plot is 
that the amplitude of the K-band outburst was considerably less (by 2.2 magnitudes) than 
in the optical. 

From Figs. UJ, [8] and and previously published data, we can conclude that: 

• The total duration of the eruptive event was between 847 and 932 days (~28 and 31 
months or 2.3 to 2.6 years). This value assumes the outburst began between MJD 2850 
and 2935 and ended around MJD 3782. 
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• The pre-outburst (SDSS, MJD 1135) and post-outburst (Gemini, MJD 3782) bright- 
ness of V1647 Ori are very similar. 

• The rise from the pre-outburst brightness to peak brightness took between 145 and 
230 days (due to the uncertainty in the date of the start of the eruption, see above). 

• The decline from peak brightness to pre-outburst brightness had two well-defined re- 
gions, a shallower phase from the late summer of 2004 to the fall of 2005 (~410 days) 
when V1647 Ori faded by ~1 magnitude, and a steeper phase from the fall of 2005 to 
its pre-outburst brightness in early 2006 (~130 days) when it faded by an additional 
~3.5 magnitudes. 

• Short-term variability was present during the shallow decline phase. This variability 
amounted to a peak-to-peak amplitude of ~0.7 magnitudes. 

• The 2003 eruptive event may have started from a somewhat elevated state. In Novem- 
ber 1998 (MJD 1135) the SDSS (r' and i') photometry lies between 1 and 2 magnitude 
fainter than the brightness when the major eruption phase commenced in late 2003. 
In addition, our last r' observation from February 2006 indicates that the star had 
returned to this fainter level after the outburst had ceased. This behavior is shown 
graphically in Fig. [TU] where all of the above optical i' and lc data are plotted together 
with an interpolated i' magnitude from 2006 February]^ We note also that the two 
lc observations of Briceno et al. (2004) from 1999 suggest that the source possessed 
significant variability when faint with photometry from January 1999 (1^=18. 44±0. 11 
transforms to i'=19.31) and December 1999 (I c =20.08±0.4 transforms to i'=20.95) 
exhibiting a statistically significant difference of ~1.6 magnitudes. This may suggest 
that V1647 Ori underwent a period of significant instability prior to the main eruption 
phase. 

• The optical (r') and near-IR (K) light-curves differ significantly in their amplitude with 
the optical being around 2.2 magnitudes larger. 

• Some 12 months after the decline to pre-outburst brightness, V1647 Ori remained 
optically very faint and close to its pre-outburst brightness. 



4 using the 2006 February r' magnitude and the r'-i' color from our 2006 January observation. 
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3.2. Optical Spectroscopic Features and their Evolution 

In this section, we discuss our optical spectroscopic observations over the monitoring 
period but defer discussion of the strongest feature in the spectra, the Ha line, to §3.3. 
However, we do consider the Ha line when describing the HIRES echelle spectra below as an 
introduction to the line properties. We note that all radial velocities quoted below have been 
corrected to be heliocentric and that to convert to velocities in the stars frame of reference 
the reader should use the rest velocity of the L1630 molecular cloud measured directly on 
V1647 Ori by Andrews, Rothberg, & Simon (2004) to be +10 km s -1 (see their Fig. 3). 

We show a GMOS spectrum of V1647 Ori in Fig. [TT] selected as representative of the 
15 GMOS spectra we obtained between February 2004 and February 2006. These are the 
data from 2004 October 6 since it has good signal to noise and shows most of the spec- 
tral features present. The source shows a rising red continuum with strong Ha emission 
exhibiting a P Cygni profile. There are numerous weak features in the spectrum including 
the Na D absorption at 5890 and 5896 A. Table H] lists the lines found in the spectra and 
their variability over the monitoring period. We have not measured equivalent widths for 
the lines since, as we shall discuss below, such a value is affected significantly by the highly 
variable optical continuum flux observed. Instead, we have merely specified whether the 
lines are in absorption or emission, or not present. As well as permitted emission lines of Fe, 
O, and Mg, numerous forbidden emission features, i.e. [Fe II], [O I], [S II], and [Ca II], are 
observed. In Figs. O [13], [HI [161 an d C3, we show extracted wavelength ranges around the 
Na D absorption lines, and the [S II], [O I] (6300 A), Fe (6400-6550 A), and the Ca emission 
lines (7200-7400 A), respectively. 



3.2.1. Na D absorption lines 

The Na D neutral absorption line profiles are shown in Fig. O The first observation 
of the monitoring period, taken on 2004 February 14, only showed weak Na D absorption. 
This was blue-shifted by -200 km s _1 with respect to the rest wavelengths of the lines. 
Between 2004 February 14 and March 10, the Na line profiles change significantly, they 
became stronger and considerably broader. On 2004 March 10, they were blue-shifted by - 
280 km s _1 with respect to the line rest wavelengths. Also at this time two emission features, 
presumably Na D lines, are present at 5895.3 A (+270 km s _1 with respect to the bluer Na D 
line rest wavelength of 5889.9 A) and 5899.7 A (+190 km s _1 with respect to the redder 
Na D line rest wavelength of 5895.9 A). These appear to be real features since emission is 
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also seen in the lower resolution 2004 February 18 spectrum from Briceno et al. (2004)13. 
An absorption feature is seen at 5877 A and its depth appears correlated to the strength 
of the Na D lines. This absorption feature is similar in width (~160 km s _1 ) to the main 
Na I lines and could possibly be a highly blue-shifted component of Na absorption appearing 
at -640 km s _1 with respect to the shorter wavelength Na I line at 5890 A. An alternative 
explanation for the presence of this line is that it is the He I line at 5876 A in absorption. 
Perhaps this is the more reasonable interpretation and is supported by the fact that we also 
observe He I lines at 6678 and 7766 A in absorption. 

After the summer months of 2004, the first spectrum taken was on 2004 September 9 
and showed the Na D lines even stronger than on March 10. Between September 9 and 2005 
January 8 (the last observation before the 2005 summer break), the Na lines vary in intensity, 
becoming weaker and stronger in successive months. The Na lines appear strongest on 2005 
January 8 and are blue-shifted with respect to the rest wavelength by ~-150 km s" 1 . Large 
changes occur between 2005 January 8 and the next observation on 2005 August 30 where 
the Na lines have weakened significantly and appear less blue-shifted (~-100 km s _1 ). In 
the UT 2005 September 25 spectrum, the Na lines possibly appear weakly in emission and 
very close to their rest wavelengths (-40 km s _1 ). From 2005 October 13 to the end of the 
monitoring period, the Na lines have more or less disappeared. 

Highly structured Na I 5890 and 5895 A lines are not uncommon amongst eruptive 
variables, specifically FUors. Examples include FU Ori itself (Bastian & Mundt 1985; Hart- 
mann & Calvet 1995), BBW 76 (Reipurth et al. 2003), V1057 Cyg (Bastian & Mundt 1985; 
Herbig, Petrov, & Duemmler 2003), and V1515 Cyg (Bastian & Mundt 1985). It was noted 
by Bastian & Mundt (1985) that all Na D profiles appear very similar to each other. This 
also is the case if we compare these three FUors to the profile of BBW 76 from Reipurth et al. 
(2003). The time-series of BBW 76 by Reipurth et al. (2003) is particularly interesting since 
they observed multiple minima in both Na lines which were predominantly blue-shifted by 
up to -300 km s _1 . In addition, between 1985 and 1994 the lines varied significantly in both 
width and blue-shifted velocity. We note here that the generally accepted interpretation of 
such Na I line structure is that an intense stellar wind, resulting from enhanced accretion, 
forms shell-like expanding/outflow structures which are accelerated close to the star/disk 
and slow with increasing distance (Bastian & Mundt 1985; Reipurth et al. 2003). Since our 
data do not have the spectral resolution of the BBW 76 data (R~2,000 as opposed to 20,000 
for BBW 76) we cannot resolve individual shell components. However, the large line width 
and the blue-shifted nature of the absorption implies in itself that we are seeing considerable 



5 Kindly made available to us by C. Briceno. 
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outflowing material with a significant velocity gradient. If the absorption feature at 5877 A 
is indeed a high- velocity blue-shifted Na component, then within one month (2004 February 
14 to March 10) the wind dramatically increases in density (line depth) and velocity (line 
width/range) with a maximum of over -600 km s _1 . For a spherically symmetric wind, the 
fact that we do not observe a symmetric profile, one with both blue-shifted and red-shifted 
absorption, suggests that we only have an unobscured line-of-sight towards the absorbing 
material that is expanding towards us. This is likely the result of obscuration by the circum- 
stellar /accretion disk and implies that the material lies close to the star and is accelerated, 
in this location, to high velocities (see the discussion of Bastian & Mundt 1985 for more 
details). This is consistent with the inferred inclination of the axis of the nebula (McNeil's 
Nebula) to the line-of-sight of ~60° (Acosta-Pulido et al. 2007; ABR08), since we are looking 
onto the star/disk through the nebula. 



3.2.2. [S II] emission lines 

Fedele et al. (2007) and ABR08 found [S II] emission at 6717 and 6731 A in the 
spectrum of V1647 Ori in 2006 January and 2007 February, respectively. In Fig. [T3], we 
trace the evolution of this emission from 2004 February to 2007 February. We see that there 
are weak [S II] lines present as early as 2004 March 10 and they are seen in all our spectra 
up until the end of the monitoring period. The lines fall at their rest wavelength with the 
resolution we have, and vary in both intensity and ratio with time. The ratio 6717/6731 is 
greater than unity from 2004 March 10 to 2006 January 5, and then, in our last two spectra, 
drops to less than unity. From 2005 October 13 onwards, the lines appear to be double- 
peaked with a separation of the peaks being ~40 km s _1 . However, our spectral resolution is 
insufficient to really study this structure in any detail or even to be certain of its existence. 
If it is present then perhaps we are detecting red- and blue-shifted shocks from the star. 

We have not derived electron densities, n e , from the ratio of the two sulfur lines 
(6717/6731, the ratio is only weakly dependent on temperature) from all our spectra since 
the signal to noise of the line detection is, in many cases, insufficient. However, for five obser- 
vations, 2005 October 13, 2005 November 27, 2005 December 25, 2006 January 05, and 2007 
February 21, we consider the lines well enough detected. The derived ratios for 6717/6731 
are 1.3±0.1, 1.4±0.1, 1.4±0.1, 1.4±0.1, and 0.8±0.1 which give values of n e , derived using 
the IRAF nebular .temden program for a temperature of 10 4 K, of 26-240, <121, <121, 
<121, and 910-2230 cm -3 , respectively. The ranges come from the uncertainties quoted on 
the line ratios and the upper limit comes from a failure if temden to calculate densities for 
a ratio >1.4. It seems, therefore, that the forbidden sulfur lines are formed in a rarefied gas 
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which becomes denser towards the end of our monitoring period although, even then, it is 
far from the [S II] critical density of 20,000 cm -3 . 

3.2.3. [0 1] 6300 A line 

The [O I] emission line at 6300 A varied significantly in intensity over the monitoring 
period. The changes are shown in Fig. [TU We note that the wavelength of this line is 
coincident with a strong night-sky emission line although we consider sky subtraction to be 
generally very good. At the start of the observing period the line was absent and appeared 
some time between 2004 February 14 and 2004 September 9 (unfortunately the spectrum 
from 2004 March 10 had a series of bad pixels at this wavelength and it was not possible to 
determine the presence of the line). From 2004 September 9 onwards, the [O I] line appeared 
to have an extended blue wing which persisted until the line faded. The blue wing extended 
to ~-270 km s~ x from the line rest wavelength and Fig. [TBI shows its profile from the spectrum 
taken on UT 2004 November 13. The line continued to increase in intensity to 2005 January 
8 where it was at a maximum. Subsequently, the line faded but was still perhaps detectable 
in the low-resolution spectrum taken on 2007 February 21. We note that the spectrum from 
2006 February 16 had bad pixels at 6300 A which were removed by interpolation. However, 
the wings of the [O I] line are still detected. 

The presence of forbidden emission lines in the optical spectra of young CTTSs has 
been studied by several authors (e.g. Edwards et al. 1987; Hartigan, Edwards, & Ghandour 
1995) and is generally associated with low-density outflowing gas. Typically, only blueshifted 
emission wings are observed which has led to the model that the circumstellar disk occults 
redshifted emission from a latitude-dependent wind with the highest gas velocities at the 
polar region of the star. From the models of Edwards et al. (1987), the only way to obtain 
[O I] profiles like those seen in V1647 Ori, i.e. single-peaked with an extended blueshifted 
wing, is to view the wind from a low inclination angle, specifically, z<45° for a wind opening 
angle 6 of <30° (see Fig. 10 in Edwards et al. 1987). This upper limit to the inclination 
angle is somewhat at odds with that derived by Acosta-Pulido et al. (2004) and mentioned 
above. 

Also, in Fig. [TU longward of the 6300 A [O I] line, we observe an absorption feature 
that deepens and broadens. Its wavelength is 6347A and we identify it as Si II. A weaker 
Si II line is also seen at 6371 A. The 6347 A line appears first weakly in the spectrum from 
2004 March 10 and becomes stronger until 2005 January 8. After this, it is possibly present 
until 2005 September 25 and then is absent. 
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3.2.4- emission lines 

In the wavelength range 6400 to 6540 A, there are four Fe emission lines from permitted 
(three) and forbidden (one) transitions. These are [Fe II] at 6432 A, Fe II at 6457 A, Fe I 
at 6495 A, and Fe II at 6517 A. The evolution of these lines is shown in Fig. [161 The 
[Fe II] 6432 A line is seen from the start of the monitoring period, 2004 February 14, to 
2005 November 27 with peak emission occurring on 2005 January 8. The Fe II 6457 A line is 
relatively weak in all spectra from 2004 September 9 to 2005 November 27 and peak emission 
also occurs on 2005 January 8. The Fe I 6495 A line is seen from 2004 March 10 to 2005 
October 13 and again, the emission peaks on 2005 January 8. Finally, the Fe II 6517 A line 
is visible from 2004 February 14 to 2005 November 19 and it is strongest on 2005 January 
8. The Fe I 6495 A line is consistently broader than either the Fe II or [Fe II] lines and 
it has a full-width half maximum (FWHM) of 250 km s _1 compared to 150 km s -1 for the 
permitted and forbidden ionized Fe lines. For comparison, the [Ca II] line at 7292 A has a 
FWHM of 150 km s _1 while a CuAr arc lamp line located near the [Ca II] line has a FWHM 
of 130 km s _1 . 

In Fig. [T7] we also identify an additional [Fe II] emission line at 7388 A. Its behavior 
is similar to that of the [Fe II] 6517 A line although it disappears slightly earlier (2005 
September 25 as opposed to 2005 November 27). It also peaks in intensity on 2005 January 
8. 

The behavior of the [Fe II] line at 6432 A seems consistent with that of the [S II] lines 
described above and suggests a common excitation mechanism such as shock-excitation. The 
permitted Fe lines are generally considered as indicators of chromospheric activity in T Tauri 
stars (Beristain, Edwards, & Kwan 1998) and therefore should show similar behavior to the 
[Ca II] lines described below. 

3.2.5. [Ca II] emission lines 

We have identified two [Ca II] emission lines at 7292 and 7324 A. The variations in 
these lines with time are plotted in Fig. [TT1 Our 2004 February 14 and March 10 spectra 
did not include these lines, however, they are present in the spectra from 2004 September 
9 to 2007 February 21. The [Ca II] lines are well detected from 2004 September 9 to 2005 
August 30 and are weakly detected in the 2005 September 25 to 2005 October 13 spectra. 
They are absent in the spectra taken after the latter date. The [Ca II] lines are unresolved 
with respect to the line width defined by the CuAr arc lamp lines (see above). 
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3. 2. 6. The Keck II HIRES spectrum 

The HIRES echelle spectra of VI 647 Ori shows a number of both absorption and emis- 
sion features with high-resolution providing information on their structure. Below, we con- 
sider each spectral feature and discuss their appearance. 

The Ha line at high-spectral resolution is particularly interesting. It has been previously 
observed at low-resolution to exhibit a P Cygni profile with a strong blue-shifted absorption 
and very broad emission (e.g. Reipurth & Aspin 2004; Walter et al. 2004). Our HIRES 
spectra are the highest spectral resolution data (R~46,000) obtained on the source to date. 
Fig. [TBI shows the profile of Ha from 2004 September 24, some 10 months after the eruption 
occurred. At this time the source had a magnitude of r'~18. The emission feature is asym- 
metric with an extended red wing reaching around +360 km s _1 . The blue-shifted absorption 
seems to begin close to the line rest wavelength and continues out to ~-360 km s _1 , the 
slight peak at this velocity appears to be where the blue wing of the emission line reaches 
the continuum level. This would make the emission symmetric about the line rest wave- 
length with a full-width zero intensity (FWZI) of around 720 km s _1 . The emission red 
wing appears to exponentially decay to the continuum. It is difficult to say where the line 
emission peak flux is located due to the overlying absorption. The absorption component 
appears to start close to the line rest wavelength and creates the very asymmetric profile 
seen in the blue wing. This is probably formed by a combination of two absorption features. 
The blue-shifted absorption goes below the local continuum level and is flat-bottomed from 
about -100 to -330 km s^ 1 . In Fig. [TBI we show the Ha profile together with a best-fit profile 
(top) fitting the red wing of the emission only. Also shown (bottom) is the residual profile 
showing the difference between the observed and best-fit profile. The profile of the fitted 
feature is Lorentzian (neither a Gaussian or Voigt profile fit the red wing well). We note that 
the profile of an unresolved arc line (ThAr) is best-fit with a Gaussian of FWHM~6 km s _1 . 
The residual profile shown in Fig. [18] exhibits extensive absorption superimposed on the Ha 
emission and clearly shows two distinct components, one at the Ha rest wavelength (the 
dot-dashed line) with a FWHM of ~50 km s" 1 , and the other extending from the continuum 
at -70 km s _1 out to -350 km s^ 1 . The more highly blue-shifted component has a more 
complex profile with the absorption being deepest at around -100 km s -1 then decreasing in 
an asymptotic manner to -290 km s _1 where the decrease becomes linear to the level of the 
continuum. 

We detect the K I 7665 and 7699 A lines which also possesses P Cygni profiles. The 
shorter wavelength line is located in the middle of the atmospheric O2 absorption band, 
however, the P Cygni structure is still evident and consistent with the longer wavelength 
line. The P Cygni absorption component has a velocity offset from the rest wavelength of 
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-160 km s 1 . The 2004 October low- resolution GMOS spectrum shows these K I lines which 
also display P Cygni profiles. 

The Na D lines are present in absorption only and are broad, relatively symmetric, and 
blue-shifted by about -145 km s -1 . These are consistent with our low-resolution GMOS 
spectrum from our nearest date, UT 2004 October 6. The signal to noise in the HIRES 
spectrum at the wavelength of the Na D lines in insufficient (S/N pixel -1 ~5.5) to say if a 
narrow interstellar Na D absorption component is present. 

The 6347 and 6371 A Si II lines are seen in absorption and both are shallow and broad 
(FWHM~-130 km s -1 ). These lines are also seen in our low-resolution GMOS spectra from 
2004 October 6. 

The Ca II line at 8498 and 8662 A lines are strongly in emission (their equivalent width 
are 10 and 8 A, respectively). The middle of the Ca II triple lines, at 8542 A, is unfortunately 
not within our observed spectral range. The longer wavelength line has a pronounced blue- 
wing absorption feature extending to about -200 km s -1 from the line center. The profile of 
this line is shown in Fig. [JJO The line is triangular in profile as is the 8498 A line. 

Strong line absorption at 7773 A corresponds to the O I triplet lines. This line is 
at the end of an echelle order and so no further information is obtainable. In our low- 
resolution GMOS spectra, the O I triplet is highly blended but is seen in absorption from 
2004 September 3 to 2005 October 13. After this the lines are either not present or hidden 
in the noise. 

We detect two very weak Mg II absorption lines at 7877 and 7897 A. In our low-resolution 
spectrum from 2005 October 6, these lines are seen weakly in absorption confirming their 
detection in the HIRES data. 

It seems that the 6300 and 6363 A [O I] lines are very weakly in emission. This is 
also seen in our 2005 October 6 low-resolution spectrum confirming their identification in 
the HIRES data. The signal to noise of the HIRES detection is low but the lines appear 
real. Also they are relatively narrow with respect to other emission lines with a FWHM of 
~40 km s -1 . 

A number of broad, weak emission lines of Fe I (6 lines) and Ti I (2 lines) are present. 
The FWHM of both the Fe I and Ti I line is ~100 km s _1 . There is a hint that the lines are 
double-peaked but, with the signal to noise present, it is difficult to say if this is correct or 
if they are broadened single lines or possibly overlapping pairs. 

From the above detailed description of the HIRES spectra, we can conclude that the 
features present agree well with those seen in the 2004 October 6 low-resolution GMOS 
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spectrum. The HIRES spectrum also confirms that some of the weaker lines seen in the 
GMOS spectrum (with low signal to noise) are real. 

3.3. The variability of the Ha line 

In Fig. [20], we present the Ha profile of V1647 Ori over the period February 2004 to 
February 2007. These 15 observations show considerable variation in both emission and 
absorption characteristics as the eruption proceeded and during the decline phase. The 
bottom right tile of the plot shows the instrumental profile (CuAr arc lamp lines) close to 
the wavelength of Ha. Fig. I2T1 shows an expanded view of the Ha profiles with an x-axis of 
velocity offset in km s" 1 from the rest wavelength of Ha (6562.8 A). 

Qualitatively, the Ha profile varies significantly even between consecutive observations. 
Ha typically exhibits P Cygni type structure with an emission component and a blue-shifted 
absorption feature. Throughout the outburst period, the emission component is extremely 
broad with a maximum FWZI of close to 1500 km s _1 , however, it decreases in width as the 
eruption progresses. The width of the blue-shifted absorption is over 600 km s -1 at the start 
of the monitoring period and also becomes smaller with time. Additionally, the absorption 
shifts towards the rest wavelength of Ha as well as decreasing in strength. Throughout 
the first 18 months of the monitoring period, the absorption component absorbs below the 
continuum (see panels 2004FEB14 to 2005AUG30 in Figs. M and EJ. After 2005AUG30, 
the wind absorption appears to be not strong enough to absorb below the continuum level 
and merely eats into the Ha emission. Finally, we note that the Ha emission and absorp- 
tion components are not of sufficient signal-to-noise in our last three datasets (2005DEC25, 
2006JAN05, and 2006FEB16) to reliably determine the extent of the absorption on the 
emission component. 

We quantify the Ha profile structure and variations in both Tableland Fig. [221 In the 
latter, we show several measures of Ha profile structure, specifically the equivalent width 
of the emission component (Wa, top-right), the full-width zero intensity (I max , FWZI) and 
the full- width at I m(KC /40 (2.5%) intensity (FW2.5%, middle-left), the velocity width of the 
blue-shifted absorption feature (AV, middle- right), the velocity offset of the wavelength 
of deepest or 'characteristic' absorption {V c h ar , bottom-left), and the wavelength of peak 
(i.e. Imax) Ha emission (X pea k bottom- right). We show the FWZI and FW2.5% variations 
since it is interesting to both compare the behavior of these measures of profile width and 
directly relate the latter to values seen in other stars, specifically from Reipurth, Pedrosa, & 
Lago (1996, henceforth RPL96). RPL96 pointed out that the FW2.5% is considerably less 
sensitive to noise than FWZI. 
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From Fig. [22] we see that as the eruption proceeded from early 2004 to early 2006 
(detailed in the r' band brightness variations shown in the top- left panel), the behavior of 
the Ha profile measures was rather complex. Clearly, detailed numerical modeling of the 
profile changes is required to characterize and quantify the physical and geometric properties 
of the accretion and outflow during the eruption (see for example Kurosawa, Harries & 
Symmington 2006). However, we consider that the above five empirically defined quantities 
may give us some useful insight into the nature of the variability observed. 

3.3.1. Ha equivalent width in the 'High Plateau 7 phase 

Our measurements of the variation of the Ha emission vs. time are shown as 

the filled circles in the top-right panel of Fig. [221 This demonstrates that from soon after 
the outburst began (2004 February) through to the start of the major decline phase (2005 
September), was approximately constant, with perhaps a slight decline, and had a 

mean value around -30 A. This time period was referred to as the 'high plateau' phase of 
the outburst decline by Acosta-Pulido et al. (2007) and we herein adopt this descriptive 
designation. Over the high plateau period, the brightness of V1647 Ori showed a decline of 
about ~1.2 magnitudes (see the top-left panel in Fig. [221) . This would result in an increase 
in Wa for constant Ha emission flux. Our average value of W# Q ~-30 A is consistent with 
the values presented in both Walter et al. (2004) and Ojha et al. (2006) whose data are also 
shown in Fig. [22] (as open squares and open circles, respectively). There are clearly short- 
term fluctuations in Who in all three datasets. These variations were considered possibly 
periodic in nature by Walter et al. (2004) and Ojha et al. (2006). Acosta-Pulido et al. (2007) 
detected short-term optical photometric variability in their light curves and, from a Fourier 
periodicity search, derived a 56 day period. They commented that one possible interpretation 
of such a periodicity would be variable obscuration by circumstellar disk material orbiting 
the young star. 

3.3.2. Timescales of Ha variability 

If such a periodicity is present in broad-band optical photometry or W# Q , or both, then 
it is clearly an important phenomenon for investigating the physics occurring during the 
V1647 Ori outburst. In our spectroscopic data, we perhaps see some temporal structure 
that could be interpreted as periodic in nature. For example, in the period September 2004 
to January 2005, the variations appear reasonably well phased with the data of Walter et 
al. (open squares). However, with so few observations (five) and their temporal spacing, 
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we cannot confirm that such periodic variability exists. What we can do is characterize the 
shortest timescale upon which variability is found. To do this, we consider the minimum 
separation between consecutive observations (combining all data points) which show statis- 
tically significant variability and we adopt a >10% change in Wj/„ as indicating statistical 
significance. There exists in the data only one date (MJD 3378) when observations were 
taken less than 1 day of each other (i.e. on the same night). Upon this date, the Wa values 
are consistent to better than 10%. Next, there are three occurrences of consecutive observa- 
tions separated by 1 day (i.e. taken on consecutive nights) and in all three cases variability 
was not detected at the >10% level. The shortest time interval over which a >10% change 
in W Ha is detected is 3 days. Specifically, between MJD 3066 and 3069, W Ha changed by 
+6 A (from ~-40 A), a change of ~15%. We conclude that, from the limited data available, 
variability in W# Q is present in VI 647 Ori on timescales as short as 3 days. 



3. 3. 3. Ha equivalent width in the decline phase 

After the period of relatively constant Wh q (February 2004 to September 2005), Wn a 
begins to change dramatically as the major photometric decline phase starts. By the time the 
source has faded to close to its pre-outburst brightness (February 2006), Wj/« has increased 
to about -100 A. The most obvious cause for such a dramatic change is the decrease in 
continuum emission as the source faded. From September 2005 to February 2006 the observed 
r' band magnitude decreased by ~4 magnitudes or a factor 40. If the Ha flux remained 
constant during this period then we would expect Wn a to have increased by a factor 40 
from ~-30 A to -1200 A. Since the increase in W Ha was measured to be only a factor 
~3.3, we can conclude that the Ha emission flux must have significantly decreased from 
September 2005 to February 2006. We estimate this decline to be about a factor ~12. If the 
Ha emission is proportional to the mass accretion rate then this decline would correspond 
to an order of magnitude decrease in accretion from September 2005 to February 2006. 
Using the relationship between near-IR H I emission line flux and accretion rate defined by 
Muzerolle et al. (1998a), Acosta-Pulido et al. (2007) estimated that the mass accretion rate 
declined from ~5xl0~ 6 M yr" 1 in March 2004 to ~5xl0~ 7 M yr" 1 in May 2006. This 
factor 10 reduction in accretion rate is consistent with the 12 x reduction in Ha emission 
flux we observe between similar dates. 
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3.3.4- Full Width Zero and 2.5% Intensity 

The FWZI of Ha emission can be considered a measure of twice the maximum velocity 
Ha emitting gas attains in the process of line creation. However, we note that this quantity 
is not just affected by bulk gas motion, since line broadening can occur. Muzerolle, Calvet, 
& Hartmann (1998) suggested that Stark broadening could be an important mechanism in 
producing very large line widths. Muzerolle et al. (2001) presented model Ha line profiles 
under the influence of Stark broadening by considering a realistic range of physical parame- 
ters related to the accretion process, magnetosphere size, inclination angle, and stellar mass 
and temperature. A comparison of Ha line profiles both with and without Stark broadening 
(and continuum opacity) showed that this effect can be a significant contributory factor with 
an increase of up to 2.5 x in FWZI. Their consideration of such broadening resulted in line 
profiles that better matched those observed in T Tauri stars. Stark broadening requires 
optically thick conditions and very high densities (n# <10 12 cm -3 ), both of which are found 
in post-shock regions of hot accretion shocks (Hartigan et al. 1991; Valenti, Basri, & Johns 
1993). 

As mentioned above, one problem with studying FWZI is signal to noise (henceforth 
S:N). Poor S:N can mask the true wavelength at which line emission fades to the continuum 
level. In the case of V1647 Ori, this is particularly important in the later spectra when 
the source had faded below r'=20. Table [5] shows the S:N at both l max and in the adjacent 
continuum for all spectra. FWZI is likely a good measure of velocity width from February 
2004 to the summer of 2005 where the S:N in the continuum is always >20. However, after 
this time the continuum S:N drops to <5. Alternatively, FW2.5% of Ha emission has been 
considered (by RPL96) to be a measure of line width less sensitive to S:N since it does not 
require a precise definition of continuum level. This is due to FW2.5% being of a larger signal 
level and hence insensitive to continuum S:N issues. RPL96 observed a significant number 
(63) of young low- and intermediate-mass stars at high spectral resolution (R~50000) and 
presented Ha line profiles for each. The profiles were divided into four classes (Type I-IV) 
depending on the asymmetry of the line profile, and the location and depth of any associated 
absorption. The reader is referred to RPL96 for the detailed definitions of the classes. 

The changes occurring in FWZI and FW2.5% with time (middle-left panel of Fig. [221) 
mirror each other reasonable well from February 2004 to the middle of 2005. This suggests 
that our measurement of FWZI is most likely accurate. In the above time period, FW2.5% 
is ~1.4x smaller than FWZI. Subsequently, however, the FWZI and FW2.5% values become 
much more similar suggesting that our values of FWZI are likely underestimates due to poor 
S:N and that the trend seen in FWZI is not purely intrinsic to the star. If we measure the 
FWZI and FW2.5% of a pure gaussian profile we find that FW2.5% is a factor 1.515 smaller 
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than FWZI. Hence, the red wing of Ha in V1647 Ori between February 2004 and mid-2005 
have FWZI and FW2.5% values of close to the gaussian ratio value. We have therefore 
multiplied the FW2.5% values from August 2005 onwards by a factor 1.515 and these are 
displayed in Fig. [22] (middle-left panel) as the FWZI values. 

FW2.5 % decreases during 2004 from a maximum close to 1000 km s _1 (FWZI~1400 km s _1 ) 
to ~750 km s _1 (FWZI~1100 km s -1 ). After this time, the decrease in FW2.5% ceases and 
from 2005 onwards it remains at ~750 km s _1 (scaled FWZI~1150 km s _1 ). Hence, through 
the major decline phase of V1647 Ori until it reaches its pre-outburst brightness, FW2.5% re- 
mains approximately constant. We note that the values encountered in V1647 Ori throughout 
the eruption are considerably larger (by a factor ~3 x ) than those found in CTTSs suggesting 
that accretion activity has not declined completely to a quiescent level. 

The range of values of FW2.5% seen in V1647 Ori is 750-1500 km s _1 . Fig. 9 of 
RPL96 shows that of their 63 sources, the peak in the distribution of (red wing) FW2.5% 
values lies at ~600 km s _1 with only three stars showing values greater than 1000 km s _1 . 
Hence, even when V1647 Ori had faded to its pre-outburst brightness, its FW2.5% value 
remained larger than the RPL96 distribution peak. The three sources that RPL96 found to 
have FW2.5%>1000 km s -1 were R Mon, FW Cha, and AS 353A. These sources were all 
classified as III-E0 in the RPL96 scheme possessing blue-shifted absorption. In total, 15 of 
their 63 sources showed type III-B profiles. The Ha profiles of V1647 Ori suggests that its 
RPL96 classification evolved with time from i) IV- E@ in February/March 2004, to ii) III-B in 
September/October 2004, to iii) IV-B in November/December 2004, to iv) III-B in January 
2005 through to at least October 2005 and possibly beyond to February 2007 (although the 
S:N of the spectra are too poor to be definite). 



3.3.5. Absorption width and velocity offset 

The width of the blue-shifted absorption feature (herein termed AV) and the velocity 
of the center of the blue-shifted absorption (herein termed v c har) gives us information on the 
velocity field of the outflowing absorbing gas. Specifically, we consider that we can associate 
AV with the velocity dispersion in the ejected material, and v c h ar with its 'characteristic' 
velocity. 



6 Double-peaked Ha emission with the peaks being of unequal height and the secondary peak being <50% 
of the primary peak 

7 P Cygni profile with absorption at sufficiently high velocity absorption to absorb continuum flux beyond 
the wing of the Ha emission. 
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The temporal variability of AV (middle-right panel of Fig. [221) is considerably different 
from that of the Ha emission component. A very rapid decline in absorption width is seen 
between our first two observations taken in February 2004 (AV~600 km s _1 ) and March 
2004 (AV~300 km s _1 ). After this time, AV is approximately constant to the end of our 
winter 2004-2005 monitoring season. At the start of the fall 2005 observing period, AV has 
decreased to ~100 km s -1 and there is some suggestion that this decrease was related to the 
main photometric decline. During the fall/winter months of 2005, AV remains approximately 
constant although the S:N of the spectra is poor. The final observation point, February 2007, 
also had very low S:N but the indication is that the value of AV had changed little. 

The velocity offset of the deepest absorption, v c har (bottom- left panel, dashed line), mir- 
rors AV showing a similar large change between the first two observations (from ~500 km s _1 
to ~300 km s _1 ). This is followed by a slower decline (to ~200 km s -1 ) between the spring 
and fall of 2004. In the fall of 2005 v c har had again decreased, this time by a factor ~2 (to 
<100 km s _1 ). By early 2006, when V1647 Ori had returned to its pre-outburst brightness, 
the S:N of the spectra and the weakness of the absorption made a determination of v c har 
difficult and all we can say is that it was <100 km s _1 . 



3.3.6. Peak emission wavelength 

The wavelength of peak Ha emission is a rather complex quantity since it is dependent 
not only on the motion of the emitting material, but also on the characteristics of any 
associated absorption. For example, if a blue-shifted wind absorbs the blue wing of the 
emission line then the wavelength of peak emission can shift to the red. Alternatively, if the 
blue-shifted absorption is weak and we measure a wavelength shift in peak emission, then 
we could interpret this as evidence for bulk motion of the emitting gas. Multiple absorption 
components, as we saw in Fig. [TBI also complicates the interpretation. In light of the above 
discussion, we refrain from any ad hoc interpretation of the changes in X pea k and merely 
relate the variations observed. 

The final panel of Fig. [22] shows the wavelength of peak Ha emission (\ pea k, bottom- 
right panel) over the monitoring period. Here we see a peak-to-peak variability of ~3 A 
(~140 km s _1 ) from February 2004 to February 2007. Soon after outburst ( February/March 
2004), the peak in emission was measured to be slightly to the red of our assumed Ha rest 
wavelength (6562.8 A). By the fall of 2004, X pea k had shifted to the blue by about 1 A 
(~45 km s _1 ). During the fall/winter of 2004, the peak had shifted back to its original 
wavelength (of February/March 2004). When our monitoring commenced again in the fall 
of 2005, \p e ak was at about the same wavelength as in spring 2005 but soon after it shifted 
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briefly to the blue and then rapidly to the red and by the spring of 2006 it was located 
redward of the rest wavelength by ~1.7 A (~80 km s _1 ). This trend continued until our last 
observation, in February 2007, when \ pea k was at +2.5 A (~120 km s _1 ). 

3.4. Ha Profile Modeling 

It is generally agreed that the presence of complex Ha emission structure with blue- 
shifted absorption in the spectrum of a young star is the result of the combination of emission 
from magnetospheric accretion and absorption from strong, dense stellar winds. The recent 
papers by Muzerolle, Calvet & Hartmann (1998), Alencar et al. (2005), Kurosawa, Harries, 
& Symington (2006), and Bouvier et al. (2007) discuss models involving such processes and 
their fit to observations of low-mass stars. It is clearly important to consider such numerical 
simulations to available datasets to allow the derivation of physical parameters such as 
outflow and accretion rates, and test the validity of the specific models. However, due to the 
complexity of such analyses in regions of high accretion (Kurosawa, private communication), 
what follows will concentrate on a qualitative understanding of the changes observed. 

4. Physical Components of the Outburst 

The eruptive event that commenced between 2003 October 23 and 2003 November 15 
(Briceho et al. 2004) resulted in the dramatic brightening of V1647 Ori and the appearance 
of what is now known as McNeil's Nebula. The event was apparently extremely violent 
since about three months after the outburst occurred (when our first observations were 
taken) the star had optically brightened by a factor ~100 and produced a P Cygni Ha 
profile exhibiting both a very high velocity (~600 km s^ 1 ), blue-shifted, absorption trough 
(indicative of the presence of a dense wind) and strong emission with an enormous velocity 
width (i.e. FWZI~1500 km s _1 - indicative of magnetospheric accretion and perhaps Stark 
broadening). How these Ha spectral features evolved over the first three months of the 
eruption is unknown, however, they are likely to be at least similar in nature to those 
observed at three months old if not even more extreme. Below, we summarize what we can 
deduce from our results regarding the properties and evolution of the wind and accretion 
phenomena. 
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4.1. The Fast Wind 

At the three month mark of the eruption, a strong wind existed producing two minima 
in the blue-shifted absorption (see Fig. [21]). Such a double-trough absorption profile was 
evident at times, for example, in the spectra of V1057 Cyg, a classical FUor (see Bastian 
& Mundt 1985). However, the two components of the wind disappeared between our first 
and second observations on UT 2004 February 14 and UT 2004 March 10, respectively. 
The optical spectrum of Ojha et al. (2005) from UT 2004 February 22 did not show such 
features either. The two minima in the Ha absorption on UT 2004 February 14, were at 
velocities of v= -270, and -530 km s _1 with respect to the rest wavelength of Ha (6562.8 /xm). 
On UT 2004 March 10, the single absorption minimum was at v= -350 km s _1 . A double- 
peaked absorption could be associated with multiple shell-like wind components ejected with 
different velocities or with the geometry and motion of the emission region. A continuing 
slowing of the wind is indicated by the v c har values plotted in Fig.[22lwhere the velocity of the 
absorption minimum declined from a maximum of -530 km s~ l to a minimum of -150 km s" 1 . 
The width of the absorption component, AV, also declined over this period (from 560 to 
100 km s _1 ) indicating a reduction in velocity dispersion in the wind and suggesting that the 
wind decelerated as it expanded away from V1647 Ori. Such behavior is perhaps consistent 
with the stochastic wind models of Grinin & Mitskevich (1992) and Mitskevich, Natta, & 
Grinin (1993). In these models the wind is non-isotropic and clumpy and decelerates at large 
distances from the young star. 

4.2. The Emission Region 

During the period of rapid decline of v c har, he. between 2004 February 14 and March 
10, the FWZI of the Ha emission remained relatively constant at ~1500 km s^ 1 . As we saw 
above, such a value is somewhat larger than typically found in accreting classical T Tauri 
stars (from the work of RPL96). Additional examples of this are RW Aur, which exhibited a 
FWZI of ~1000 km s _1 in the spectra of Alencar et al. (2005), and some sources studied by 
Muzerolle, Hartmann, & Calvet (1998) and Alencar & Basri (2000) where a typical FWZI 
was found to be in the range 600 to 1000 km s _1 . The FWZI did not show significant 
decline until sometime between our observation in March 2004 and our next observation 
in September 2004. Over this period the FWZI had reduced by around 500 km s _1 (to 
~1000 km s _1 ). During the period of the major photometric decline, from August 2005 to 
February 2006, the FWZI remained in the range 1000-1400 km s _1 , i.e. the dramatic change 
in optical brightness of V1647 Ori did not appear to be strongly correlated with the FWZI 
of the Ha emission. This argues that the optical fading and the Ha emission are not the 
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result of the same physical processes. 

As we related above, contributions to the FWZI of the Ha emission probably arise in 
both bulk gas motion and Stark broadening. The region of Ha emission therefore much be 
extremely dense and suggests that it is the post-shock region of the accretion shock. The 
change in optical colors as the major photometric decline phase began (they became redder) 
suggests that there was a reduction in blue continuum from the accretion shock region. 
Briceno et al. (2004) observed the spectrum of V1647 Ori using reflected light from McNeil's 
Nebula and found it to correspond to that from an early B-type source i.e. the hot accretion 
shock. Correspondingly, the reduction in blue continuum likely resulted from of a decline in 
post-shock temperature and density, probably due to a slowing of the accretion rate. The 
decline in accretion rate was estimated, using near-IR H lines, to be an order of magnitude 
between March 2004 and May 2006 (Acosta-Pulido et al. (2007). 



5. Outburst Timeline 

We now attempt to piece together the timeline of events over the outburst period and 
relate the photometric and spectroscopic variations discussed above. 

1. Prior to the outburst, V1647 Ori appeared to be in an unstable state with variations in 
optical brightness of up to ~2 magnitudes. It appeared that the main outburst event 
was initiated from an photometrically elevated state existing for at least 5 years prior 
to the main eruption i.e. from the first Briceno et al. (2004) and McGehee et al. (2004, 
SDSS) observations. 

2. The outburst began between late- July and late-October 2003. Unfortunately, there are 
only a handful of photometric measurements for the first three+ months of the event. 
These were obtained serendipitously during an Orion OBI photometric monitoring 
program (Briceno et al. 2001). The rise time to peak brightness was short, 120 days, 
with a change in optical (I-band) brightness of 3.6 magnitudesjfl 

3. At the time of our first observations, the first taken on the source after its discovery by 
McNeil (2004), the star had brightened from its pre-outburst level by over 5 magnitudes 
in the optical and over 2 magnitudes at 2 /im. It had optical spectral features indicative 
of a massive accretion burst producing a fast, dense stellar wind creating high- velocity 



8 V1647 Ori had already brightened by ~1 magnitude prior to this but it is unknown when this occurred 
due to lack of observations during this period (Briceno et al. 2004). 
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blue-shifted P Cygni-type absorption on strong Ha emission. In the NIR, VI 647 Ori 
possessed strong emission in Br7, Na I, and the v = 2 — CO overtone bandheads, all 
indications were that the emission region was hot, extremely dense, and in an excited 
kinematic state. The optical brightening is attributed to a combination of the star 
being revealed after dust in the immediate circumstellar environment of the object was 
sublimated (Aspin et al. 2009) together with a luminosity increase due to the addition 
of significant hot accretion flux. The brightening in the NIR is perhaps more likely the 
effect of irradiation and subsequent heating of the inner regions of the circumstellar 
disk, beyond the sublimation radius, by UV flux from the accretion flow. However, it 
is clear that the heating did not propagate through to outer regions of the disk since 
the sub-mm/mm flux remained unchanged. 

4. Over the next 10 months, V1647 Ori remained within half a magnitude of its outburst 
optical/NIR brightness although there was seemingly stochastic variability at the level 
of around 1 magnitude peak-to-peak. Statistically significant brightness fluctuations 
were seen on periods as short as a few days. During this so-called "high-plateau" phase, 
the optical spectrum of V1647 Ori showed strong P Cygni Ha emission on a relatively 
stable continuum. This suggests that the accretion process continued unabated over 
this period with only minor variations occurring. This is supported by the reasonably 
constant FWZI of Ha (~1500 km s^ 1 ), indicating that the emission was taking place 
in very dense gas, and the derived accretion rate from NIR line emission. During 
this same period, the wind quickly slowed (from v c har ~600 km s _1 to 200 km s _1 ). 
Also, the velocity dispersion in the wind mirrored v c h ar in that it quickly declined (from 
AV~600 km s -1 to 200 km s -1 ). During the rapid decline in wind velocity and velocity 
dispersion, weak shock-excited emission lines of [S II] appeared in the optical spectrum. 
These remained approximately constant to the end of the monitoring period. During 
the summer months of 2004, an [O I] emission line appeared and strengthened for 
about a year. This may be the result of the formation of a partly obscured bipolar- 
like outflow and its formation approximately corresponded to the period of the high- 
plateau. Subsequently, the [O I] line declined in intensity which corresponded to the 
rapid decline in brightness of V1647 Ori. 

5. At the start of the rapid decline phase of the optical flux (in the fall of 2005, some 
6 months after our last plateau-phase observation), V1647 Ori had already shown a 
significant reduction in Ha FWZI (~1500 km s _1 to ~800 km s" 1 ), wind velocity 
(v c fear~300 km s _1 to ~200 km s" 1 ), wind velocity dispersion (AV~250 km s _1 to 
~150 km s" 1 ), and accretion rate (a factor 10). This suggests that the density of the 
region of Ha emission declined and the wind slowed. Two of the three aforementioned 
parameters (FWZI and AV) remained approximately constant to the end of our mon- 
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itoring period (in 2006 February) when the source had faded by ~4 magnitudes and 
returned to its pre-outburst optical brightness. The only significant change was in 
wind velocity which continued to decline (<5v~250 km s _1 to ~100 km s _1 ) into early 
2006 implying that the wind was still decelerating. 

6. In the period after February 2006, beyond the time period investigated here, we ob- 
tained one further set of observations that proved rather interesting. These were pub- 
lished in Aspin, Beck, & Reipurth (2008) and were obtained approximately one year 
after the last observation presented above. In February 2007, the optical brightness of 
V1647 Ori remained close to its February 2006 value with r'=23.3. Little evidence was 
seen for blueshifted Ha absorption suggesting that the wind had continued to weaken 
over the intervening year. The FWZI of Ha remained at ~800 km s _1 implying that 
the accretion process had perhaps settled into a quasi-steady state. The Ha line pro- 
file calculations by Muzerolle, Calvet, & Hartmann (2001) discussed above, gave values 
of FWZI of this order (see their Figs. 8, 11, and 12) and produced a good match to 
observations of some CTTSs (e.g. BP Tau with FWZI~800 km s _1 ). In addition, the 
February 2007 data showed that shock-excited [S II] emission was still present, as was 
weak [O I] emission with a weak blueshifted component. 
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Table I. V1647 Ori Gemini GMOS-N Observing Log 



UT Date 


MJD a 


Filters and/or 


Exposure Times 


Seeing 


yymmdd 




Grism Used 


(seconds) 


(r'-band) 


981116 


1134 


■ b 




0'.'9 


040214 


3049 


s' r' i' 7 ,c R831 d 


60 60 60 60 


0'.'54 


040310 


3074 


r',B600° 


10,1200 


0'.'96 


040903 


3251 


g',r',i',z' 


30,30,30,30 


0'.'53 


040909 


3257 


B600 


900 


Of.' 62 


040924 


3272 


HIRES 


3600 


? 


041006 


3284 


g',r',i',z',B600 


30,30,30,30,900 


0'.'48 


041113 


3322 


g',r',i',z',B600 


30,30,30,30,1200 


0'.'49 


041212 


3351 


g',r',i',z',B600 


30,30,30,30,1200 


Of.' 58 


050108 


3378 


g',r',i',z',B600 


30,30,30,30,1200 


0'.'78 


050830 


3612 


g',r',i',z',B600 


30,30,30,30,900 


0'.'63 


050921 


3634 


g',r\i',z' 


30,30,30,30 


Of.' 48 


050925 


3638 


B600 


900 


Of.' 65 


051013 


3656 


g',r',i',z',B600 


30,30,30,30,900 


0'.'68 


051119 


3692 


r',B600 


30,900 


Of.' 60 


051127 


3701 


B600 


30,1200 


Of.' 61 


051128 


3702 


g',r',i',z' 


60,60,60,60 


Of.' 59 


051225 


3729 


g',r',i',z',B600 


60,60,60,60,1200 


Of.' 59 


060105 


3740 


g',r',i',z',B600 


60,60,60,60,1200 


Of.' 49 


060216 


3782 


r',B600 


120,3600 


Of.' 42 


061222 f 


4090 


R 


30 


Of.' 90 


070221 


4152 


R400 


3.5 hours 


Of.' 76 


070222 


4153 


g\r',i\z' 


600 


Of.' 60 



a Modified Julian Date. 2450000+ 

b SDSS observations taken from the galactic plane Orion release data 
(Finkbeiner ct al. 2004). 

c GMOS "Sloan" filters. 

d GMOS Red 831 lines/mm grating, blaze wavelength 757nm, R~4396 
with 0'.'5 slit, simultaneous wavelength coverage 207nm, 0.034 nm/pixel. 

°GMOS Blue 600 lines/mm grating, blaze wavelength 461nm, R~1688 
with 0'.'5 slit, simultaneous wavelength coverage 276nm, 0.045 nm/pixel. 

f Taken on University of Hawaii 2.2m telescope. 



Table 2. VI 647 Qri and comparison star coordinates and photometry a 



Source 


R.A. 


Decl. 


g' 


r' 


i' 


z' 




(J2000) 


(J2000) 


(Mags) 


(Mags) 


(Mags) 


(Mags) 


V1647 Ori 


05 46 13.1 


-00 06 05 


24.74±0.76 


23.04±0.22 


20.81±0.05 


18.80±0.04 


SI 


05 46 13.0 


-00 08 15 


23.09±0.15 


21.14±0.04 


19.47±0.02 


18.37±0.03 


S2 


05 46 12.3 


-00 08 08 


23.29±0.18 


21.19±0.05 


18.74±0.02 


16.93±0.02 


S3 


05 46 11.4 


-00 07 55 


23.26±0.17 


21.39±0.05 


20.71±0.05 


19.65±0.07 


S4 


05 46 16.3 


-00 06 52 


21.93±0.06 


20.30±0.02 


18.58±0.02 


17.54±0.02 


S5 


05 46 11.6 


-00 06 28 


23.74±0.27 


21.00±0.04 


18.96±0.03 


17.22±0.02 


S6 


05 46 09.6 


-00 03 31 


23.26±0.18 


20.35±0.02 


18.33±0.02 


16.85±0.02 



a SDSS photometry taken from the galactic plane "Orion" release data (Finkbeiner et al. 
2004) dataset 'calibImage-259-r3-0534' from 1998 Nov 17 (MJD 51134). Errors are from 
quoted inverse variance measurements on the PSF fluxes. 



Table 3. V1647 Ori Optical Photometry 



MJD a g'±err(dg') r'ierr(dr') i'±err(di') z'±err(zg') g'-r' r'-i' i'-z' Qriz h 

(mags) (mags) (mags) (mags) (mags) (mags) (mags) (mags) 



1135 c 










23 


.04±0.22(6. 


.34) 


20. 


81±0 


.05(5. 


.91) 


18. 


80±0 


.04(4. 


.41) 


1 


.70 


2 


.23 


2 


.01 





,25±0. 


23 


3049 


20. 


05±0. 


15(0. 


00) 


17 


70±0.09(0. 


00) 


15 


90±0 


.08(0 


.00) 


14. 


39±0 


.06(0. 


.00) 


2 


.35 


1 


.80 


1 


.51 





31±0. 


20 


3074 










17 


,91±0.09(0. 


21) 




































3251 


20. 


06±0. 


.15(0. 


.01) 


17 


.81±0.09(0. 


11) 


15 


93±0 


.08(0. 
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a Modified Julian Date. 2450000+ 

b Rcddcning invariant colors as denned by McGehee et al. (2004). 

C SDSS photometry taken from the galactic plane Orion release data (Finkbeincr ct al. 2004) from MJD51134 

d i' photometry derived from r' and r'-i' color from MJD 3740. 

°R photometry from image taken at University of Hawaii 2.2m telescope. 
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emission, A=absorption, WE=weak, P=P Cygni profile, X=bad pixels at location of line, ?=inconclusive, -=no line evident, 0=outside observed spectral range. 
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Table 5. 


V1647 Ori 


Ha Line Profiles 








MJD a 


S:N b 


w A c 


\ d 

"peak 

(A) 


(A(km s- 1 )) 


FWZI f 


FW2.5% S 


FW10% h 


AV 








(A) 


(km s _1 ) 


(km s ) 


(km s _1 ) 


(km s _1 ) 


(km s" 1 ) 


3049 


252(34) 


-31.4+2.0 


6563.2+0.1 


+0.4(18) 


1392+30 


968+30 


534+30 


552+30 


516+20 


3074 


400(39) 


-50.9+2.0 


6563.0+0.1 


+0.2(9) 


1384+30 


1030+30 


562+30 


270+30 


294+20 


3257 


282(42) 


-24.1+2.0 


6562.3+0.1 


+0.5(23) 


958+30 


664+30 


460+30 


308+30 


230+20 


3273 


352(10) 


-18. Oil. 


6563.3+0.1 


+0.5(24) 


960+10 


592+30 


440+10 


255+20 


210±10 k 


3284 


278(36) 


-29.5+2.0 


6562.5+0.1 


-0.3(-14) 


1208+30 


822+30 


564+30 


273+30 


220+20 


3322 


261(41) 


-22.1±2.0 


6562.7+0.1 


-0.1(-5) 


1002+30 


660+30 


450+30 


332+30 


208+20 


3351 


212(29) 


-29.2+2.0 


6562.8+0.1 


0.0(0) 


1056+30 


778+30 


526+30 


248+30 


213+20 


3378 


361(49) 


-27.1±2.0 


6563.2+0.1 


+0.4(18) 


1180+40 


856+30 


556+40 


238+30 


167+20 


3612 


160(23) 


-21.2±2.0 


6563.0+0.1 


+0.2(9) 


840+40 


740+30 


476+40 


158+30 


150+20 


3638 


108(13) 


-29.2±2.0 


6562.5+0.1 


-0.3(-14) 


930+50 


714+30 


464+50 


112+30 


144+20 


3656 


141(17) 


-30.2+2.0 


6563.4+0.1 


+0.6(27) 


1190+50 


864+30 


512+50 


95+30 


142+20 


3692 


38(4) 


-37.2+2.0 


6563.6+0.1 


+0.8(37) 


970+80 


858+30 


518+80 


109+60 


127+50 


3701 


29(2) 


-64.2+2.0 


6563.4+0.1 


+0.6(27) 


688+80 


604+30 


458+80 


98+60 


108+50 


3729 


40(3) 


-64.9+2.0 


6563.7+0.1 


+0.9(41) 


876+100 


704+30 


536+100 


<100 


_ l 


3740 


39(3) 


-66.6+2.0 


6563.9+0.1 


+1.1(50) 


868+100 


686+30 


598+100 


<100 




3782 


29(1) 


-98.7+2.0 


6564.4+0.1 


+1.6(73) 


770+100 


730+30 


632+100 


<100 




4152 m 


15(1) 


-136.2±2.0 


6565.4+0.1 


+2.6(119) 


822+100 


762+30 


756+100 


<100 





a Modified Julian Date. 2450000+ 

b Signal to noise on the Ha emission peak and, in parentheses, in the adjacent continuum. 

c Equivalent width of emission component calculated from the continuum levels on the blue and red edges of the main Ha emission 
peak. The uncertainties are estimated from the spread of repeat measurements. 

d The wavelength of peak Ha emission, \ pea k 

°Offsct of peak Ha emission from 6562.8 A in A and, in parentheses, km s — 1 . 

f The Full- Width Zero Intensity of the Ha emission in km s — 1 . The velocity offset to the continuum on the red side of 6562.8 A 
was measured then doubled to obtain the FWZI value. The uncertainties are estimated from the spread of repeat measurements. 

g The Full- Width Intensity of the Ha emission in km s _1 measured to a signal of I m aa;/40 (2.5%), where l m ax is the peak line 
intensity. The width is measured on the red side of 6562.8 A and was doubled to obtain the FW2.5% value. The uncertainties are 
estimated from the spread of repeat measurements. 

h The Full- Width Intensity of the Ha emission in km s — 1 measured to a signal of 10% of the peak intensity, lmax- The width is 
measured on the red side of 6562.8 A and was doubled to obtain the FW10% value. The uncertainties are estimated from the spread 
of repeat measurements. 

'The width of the Ha absorption feature in km s _1 . The uncertainties are estimated from the spread of repeat measurements. 

^The velocity offset (in km s~ x ) of the deepest absorption from 6562.8 A, the rest wavelength of Ha. The uncertainties are estimated 
from the spread of repeat measurements. 

k Keck/HIRES spectrum. Profile not shown in Figs. 1201 and 1211 

'Emission and absorption to weak to accurately measure at signal to noise present. 
m Data from Aspin et al. (2009). 
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Fig. I.— The region of L1630 including V1647 Ori and McNeil's Nebula. This r' band image 
was taken on Gemini North using GMOS on UT 2004 February 14. V1647 Ori and the six 
calibration sequence stars, labeled SI to S6, are identified. North is at the top, East to the 
left. The scale of the image is O'.'144/pixel and the exposure time was 40 seconds in seeing 
of FWHM 0752. 



3 1 1 1 5T 1 1 1 r 

o o o 

o o o 

^ L c\2 C\2 C\2 

z' 

•- - 

i' 



00 



3 



P0 



CV2 

c\2 



r' 

1 PO 



• • • 



r 



* r 



■8 z' D „ • 



J „. ^ • 

i- % * > ^ • ■ 

\ ^ 

* 

^ |_ V 
CV2 ' * 



3000 3200 3400 3600 3800 

Modified Julian Date (2450000+) 



Fig. 2. — Optical photometry of V1647 Ori spanning the period 2004 February to 2006 
February. The data were obtained using the Gemini facility imager/spectrograph GMOS- 
N in SDSS filters g', r', i', and z'. The horizontal lines are the pre-outburst brightness of 
V1647 Ori from the SDSS Orion survey data and additionally presented in McGehee et al. 
(2004). 
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Fig. 3.— Optical photometry at g', r', i', and z' of the field star labeled S4 in Fig. [TJ 
The observations are shown as open circles with associated error bars. These values were 
derived using the SDSS photometry of all six calibration stars (see Fig. [I]) as photometric 
reference. For each passband, the associated la variation on the mean of all observations 
is shown assuming that the star S4 is not intrinsically variable. The dot-dashed lines show 
the mean value of the aforementioned data points, while the dashed lines show the magni- 
tude of S4 derived using the GMOS zeropoints (ZP(g')=27.93, ZP(r')=28.18, ZP(i')=27.90, 
ZP(z')=26.77 from Jorgensen 2009) from early 2004. We note that the difference between 
the two lines is, in all cases, smaller than the associated photometric errors. 
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Fig. 4. — Optical colors of V1647 Ori spanning the period 2004 February to 2006 February. 
The data were obtained using the Gemini facility imager /spectrograph GMOS-N in SDSS 
filters g', r', i', and z'. The dashed, dot-dashed, and dotted lines joining the observations 
are the g'-r', r'-i', and i'-z' colors, respectively. The short dashed, dot-dashed, and dotted 
horizontal lines from MJD (50000+) 3800-3900 are the pre-outburst colors of V1647 Ori 
from the SDSS Orion survey data and additionally presented in McGehee et al. (2004). 
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Fig. 5. — A color-color plot of the r'-i' and i'-z' colors of V1647 Ori shown in Fig. HI The 
solid line is the locus of stellar main sequence dwarfs taken from Finlator et al. (2000). 
Representative spectral types are shown. The dashed lines extending from the ends of this 
locus are reddening vectors for R=3.1. These are calculated using the tabular data provided 
by D. Finkbeiner (http://www.astro.princeton.edu/~dfink/sdssfilters/). The three points 
labeled (MJD 50000+) 3702, 3729, and 3740 represent data taken close to the end of the 
outburst. Their location is consistent with either increased reddening with respect to the 
outburst colors, or a change in effective temperature. The pre-outburst colors of V1647 Ori 
are labeled SDSS 1135. 
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Fig. 6. — A plot of the reddening invariant color, Q r j 2 , vs. MJD. Within the associated 
errors, the only possible trend in Q r ^ is a slight increase in value after MJD 3600. See 
Section 3.1.2 for the definition of Q r iz- 
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Fig. 7. — Optical i' band photometry of V1647 Ori using GMOS (open circles) from 2004 
February to 2006 February together with the I c photometry from Briceno et al. 2004 (filled 
circles) from 2003 November to 2004 February, and the Acosta-Pulido et al. (2007) I c pho- 
tometry (filled stars) from 2004 February to 2006 February. The Briceno et al. photometry 
has been transformed from Ic to SDSS i' magnitudes using the transformation equations 
presented in Ivezic et al. (2007). Additionally the Briceno et al. data points from MJD 
2989 onwards have had an aperture correction applied (+0.62 magnitudes) since the software 
aperture used by Briceno et al. had a radius of 4'.'1 while our aperture radius was 1". The 
horizontal black dotted line shows the November 1998 SDSS i' photometric magnitude of the 
source (i'=20.81). Prior to the time period displayed, Briceno et al. reported the I c magni- 
tude of V1647 Ori to be 18.44+0.11 (i'=19.31) in January 1999, and 20.08+0.3 (i'=20.95) 
in December 1999. 
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Fig. 8. — Optical r' band photometry of V1647 Ori using GMOS (open circles) from 2004 
February to 2006 February together with the R photometry from Ojha et al. 2006 (filled 
stars) from 2004 February to 2005 November. We have shifted the Ojha et al. data vertically 
by 1 magnitude to obtain good correspondence between the Gemini and Ojha et al. values 
throughout the overlap region. The horizontal dotted line is the November 1998 SDSS r' 
photometry (r'=23.04). 
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Fig. 9. — The optical r' band light-curve from Fig. [8] together with the near-IR K-band light- 
curve compiled from the literature. The r' band photometry is shown as open circles (this 
paper) and filled stars (Ojha et al. 2006). The K-band photometry is shown as open squares 
(from Acosta-Pulido et al. (2007) and Ojha et al. 2006). We have shifted the optical data 
vertically by 8 magnitude so that both light-curves can easily been compared. The horizontal 
dot-dashed line is the November 1998 SDSS r' photometry (r'=23.04). The horizontal dashed 
line is the October 1998 2MASS K-band photometry (K=10.27). The vertical dotted line is 
at MJD 3750 and is merely there as an aid in relating the two light-curves. 
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Fig. 10. — Optical photometry of V1647 Ori using a) GMOS i' band data (open circles) 
from 2004 February to 2006 February, b) the Briceno et al. (2004) I c band data transformed 
to i' magnitudes (filled dots) from 1999 January to 2004 February, c) the i' band SDSS data 
from McGehee et al. (2004) (open squares) from 1998 November to 2002 February, and 
d) the Acosta-Pulido et al. (2007) I c band data transformed to i' magnitudes (filled stars) 
from 2004 February to 2006 February. The last point of the GMOS data has been estimated 
from the r' band observation on that MJD using the r'-i' colors of the previous observation 
(r'-i'=2.2). 
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Fig. 11. — Optical GMOS spectrum of V1647 Ori from 2004 October 6. This is a repre- 
sentative spectrum of the 15 acquired between February 2004 and February 2006. In 2004 
October, V1647 Ori had a magnitude of r'=18.03 and was in its shallow decline phase. All 
identifiable absorption and emission lines are marked including the telluric O2 and 02-A 
bands. Note the prevalence of weak permitted and forbidden Fe emission lines. The major 
intrinsic feature is the strong Ha emission line with a blue-shifted absorption component 
forming a P Cygni type profile. Table H] shows the variation in spectral features over the 
monitoring period. 
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Fig. 12. — Variation in Na D 5890 and 5896 A absorption profiles over the monitoring period. 
The continuum has been subtracted in these spectra and they have not been smoothed so 
as to preserve the intrinsic noise structure for comparative purposes. The vertical dashed 
lines show the Na lines rest wavelengths. In 2004 February, the Na lines were barely visible 
but soon became very pronounced and blue-shifted with respect to the rest wavelength. The 
lines are highly variable in depth and shape until 2005 August (at the beginning of the 
steep decline phase) when they had faded considerably. For the remainder of the monitoring 
period (2005 September to 2007 February) the lines were barely present. 
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Fig. 13.— Variation in the shock-excited [S II] 6717 and 6731 A emission line intensity 
and ratios over the monitoring period. The vertical dashed lines show the [S II] lines rest 
wavelengths. In 2004 February, the [S II] lines were not present. In the spectrum taken on 
2004 March 10 however, they are seen faintly. The intensity and ratio of the lines varies 
from then to the end of the monitoring period. Again, the spectra have been continuum 
subtracted and not smoothed. 
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Fig. 14.— Variation in the [O I] 6300A emission line intensity over the monitoring period. 
The vertical dashed line shows the [O I] line rest wavelength. The line is not present on 2004 
February 14 but is strong in the spectrum taken on UT 2004 September 9. The spectrum 
from 2004 March 10 had considerable noise spikes at the wavelength of [O I] which have been 
removed by interpolation. The line intensity mostly increases from UT 2004 September 9 
to UT 2005 January 8 when it begins to decline. The spectrum from 2006 January 5 had 
some noise spikes (which have been removed) but the wings of the [O I] line can still be seen. 
The line is perhaps just detected in the UT 2007 February spectrum. Again, the spectra 
have been continuum subtracted and not smoothed. The [O I] emission line is asymmetric in 
shape with an extended blueshifted wing. We note the presence of an unidentified absorption 
feature at ~6350 A. 
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Fig. 15. — The [O I] 6300 A emission line profile from the UT 2004 November 13 spectrum. 
The vertical dot-dashed line shows the [O I] line rest wavelength. The vertical dotted lines 
from the bottom axis show the extent of the emission profile which corresponds to -270 to 
+100 km s" 1 . The sloping dashed line is the continuum level across the emission feature. 
The dotted line under [O I] shows an arc spectrum lamp line close in wavelength to [O I]. 
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Fig. 16.— The Fe lines in the 6400 to 6540 A region of the spectra. These lines are [Fe II] 
at 6432 A, Fe II at 6457A, Fe I at 6495 A, and Fe II at 6517 A. The vertical dot-dashed 
lines are at the line rest wavelengths. The lines vary in intensity together, getting stronger 
from UT 2004 February 14 (when they are very weak) to UT 2005 January 8. After this, 
they decline in intensity until UT 2005 November 27 when they are not detected. Again, 
the spectra have been continuum subtracted and not smoothed. 
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Fig. 17.— The [Ca II] lines in the 7200 to 7400 A region of the spectra. Specifically, at 
7261, 7292, and 7324 A. The vertical dot-dashed lines are at the line rest wavelengths. At 
the long-wavelength edge of the plot is an [Fe II] at 7388 A. Our UT 2004 February 14 and 
UT 2004 March 10 spectra did not cover this region of the spectra, however, from UT 2004 
September 9 to UT 2005 October 13, the [Ca II] lines are visible. After this they are not 
present. Again, the spectra have been continuum subtracted and not smoothed. 
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Fig. 18.— The Ha line profile from the HIRES spectrum of V1647 Ori. At the top is the 
observed spectrum overlaid with a model Gaussian profile that fits the red wing of the line. 
Note the wings of the profile extend to ±360 km s _1 from the nominal wavelength of Ha 
(dot-dashed line). For comparison, an arc lamp line lying close to the wavelength of Ha has 
a full-width half-maximum gaussian profile of 12 km s~ x . Below is the difference between the 
observed and model profile. There appears to be significant Ha absorption at zero velocity 
together with blueshifted absorption extending from -80 to -360 km s -1 . The minimum in 
the Ha absorption lies at 120 km s _1 . 
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Fig. 19.— The HIRES spectrum of V1647 Ori containing the Ca II line at 8662A. Note 
the trianglular profile of the line. Also present are two faint Fe I lines which seem to show 
a double-peaked nature. 
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Fig. 20. — Ha emission line profiles extracted from the optical spectra of V1647 Ori. The 
monitoring period is 2004 February to 2006 February and the date of observations is shown 
in the top-left corner. To better show the profile structure, the spectra have been scaled to 
a minimum and maximum flux of zero and unity, respectively. The Ha equivalent widths, 
and related information are given in Table 4. Of particular note is the prominent P Cygni 
profile in the earlier spectra. The bottom-right plot is of CuAr arc lines for comparison of 
Ha line widths and instrumental profile. 
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Fig. 21. — An expanded view of the Ha P Cygni profiles shown in Fig. [20l Here, the x-axis 
is velocity offset in km s" 1 from 6562.8 A. Note the change in emission width, absorption 
width and the trend for lower-velocity blue-shifted absorption with increasing time. 
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Fig. 22. — A plot of the changes in Ha profile structure vs. time. The units of the x-axes 
of all six panels are modified Julian Date (MJD). The top-left panel shows the V1647 Ori 
r' band photometric light curve. The horizontal dot-dashed line labeled r'po is the pre- 
outburst SDSS r' brightness. The top-right panel shows the Ha emission equivalent width 
(Wa) in Angstroms. The middle- left panel shows the Full Width Zero Intensity (FWZI), 
the Full Width 2.5% Intensity (FW2.5%), and the scaled (see text for details) FW2.5% 
Intensity of the emission component in km s _1 . The middle-right panel shows the width of 
the blue-shifted absorption component (AV) in km s _1 . The bottom- left panel shows the 
velocity offset (in km s _1 from 6562.8 A) of the deepest or 'characteristics' absorption {y char ). 
The bottom-right panel shows the wavelength of peak emission, termed \ pea k, in Angstroms. 
Velocity offsets (in km s _1 from 6562.8 A) are given inside the y-axis. 



